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We generate a second equation of state (EOS) of nuclear matter for a wide range of temperatures, 
densities, and proton fractions for use in supernovae, neutron star mergers, and black hole formation 
simulations. We employ full relativistic mean field (RMF) calculations for matter at intermediate 
density and high density, and the Virial expansion of a non-ideal gas for matter at low density. 
For this EOS we use the RMF effective interaction FSUGold, whereas our earlier EOS was based 
on the RMF effective interaction NL3. The FSUGold interaction has a lower pressure at high 
densities compared to the NL3 interaction. We calculate the resulting EOS at over 100,000 grid 
points in the temperature range T = to 80 MeV, the density range ub = 10~* to 1.6 fm~'^, and 
the proton fraction range 1^ = to 0.56. We then interpolate these data points using a suitable 
scheme to generate a thermodynamically consistent equation of state table on a finer grid. We 
discuss differences between this EOS, our NL3 based EOS, and previous EOS by Lattimer-Swesty 
and H. Shen et al. for the thermodynamic properties, composition, and neutron star structure. The 
original FSUGold interaction produces an EOS, that we call FSU1.7, that has a maximum neutron 
star mass of 1.7 solar masses. A modification in the high density EOS is introduced to increase the 
maximum neutron star mass to 2.1 solar masses and results in a slightly different EOS that we call 
FSU2.1. The EOS tables for FSU1.7 and FSU2.1 are available for download. 



PACS numbers: 21.65.Mn,26.50.+x,26.60.Kp,21.60. Jz 



I. INTRODUCTION 

The evolution of compact stellar objects provides a 
unique laboratory to study matter at extremely high den- 
sities, temperatures, and/or neutron fractions. Neutron 
star (NS) binary systems emitting X-ray bursts or super 
bursts provide a wealth of information on the NS crust, 
its composition, and transport properties HJ. Core col- 
lapse supernovae (SN) release gigantic amounts of gravi- 
tational energy through the emission of neutrinos. Their 
dynamics depends on detailed understanding of neutrino 
interactions in dense matter pi]. NS-NS and black hole- 
NS mergers produce strong gravitational waves that de- 
pend on the equation of state (EOS) of dense matter [3|. 
The EOS gives the pressure as a function of density, tem- 
perature and proton fraction. The historic detection of 
gravitational waves is anticipated soon, as the Advanced 
LIGO and VIRGO detectors become operational [ij. 

For the past two decades, there have been several at- 



tempts to construct complete EOS tables over wide range 
of temperatures, densities, and compositions. However 
most realistic supernova simulations have used only two 
EOS tables. The Lattimer-Swesty (L-S) equation of state 
0, that is based on a compressible liquid drop model 
with a Skyrme force, and the H. Sheii, Toki, Oyamatsu 
and Sumiyoshi (S-S) equation of state [1,01, that is based 
on a relativistic mean field (RMF) model with Thomas- 
Fcrmi approximation and variational method. Recently, 
Hempel and Schaffner-Bielich have built an EOS table 
by matching a nuclear statistical equilibrium model, at 
low densities, directly to uniform nuclear matter at high 
densities [8]. In addition, the S-S EOS has been extended 
to include hyperons d, • 
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In a series of papers [11| - |14| . two of the present au- 
thors presented a new complete EOS. This is based on 
an RMF model, NL3 |15|], to self-consistently calculate 
nonuniform matter at intermediate density and uniform 
matter at high density. Matter at low density is described 
with a Virial expansion for a non-ideal gas of nucleons 
[l6| and nuclei [l7(j. This Virial approach uses elastic 
scattering phase shifts and nuclear masses as input and 
includes Coulomb corrections that can be important for 
neutrino interactions [H, [l9|. The Virial expansion is 
exact in the low density limit. Altogether our RMF and 
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Virial EOS models cover the large range of temperatures, 
densities, and proton fractions necessary for astrophysi- 
cal simulations. This EOS table is available for download 

In this work, we use a different RMF effective interac- 
tion FSUGold [l^l to study matter at intermediate and 
high densities. This replaces the NL3 effective interac- 
tion used in our previous work. The low density EOS 
from the Virial expansion is unchanged. FSUGold repro- 
duces ground state properties of nuclei across the peri- 
odical table very well. The differences between NL3 and 
FSUGold have been extensively studied for example, in 
[2l| . It was observed that FSUGold is consistent with 
most known constraints, which include the universal be- 
havior of dilute Fermi gases with large scattering lengths, 
heavy-ion experiments that probe both the low- and high- 
density domain of the equation of state. NL3 may not 
satisfy these constraints as good as FSUGold. However, 
there are still large controversies in these theoretical and 
experimental constraints. Therefore it is very useful to 
have EOS produced with both effective interactions. The 
astrophysical simulations with the two EOS could poten- 
tially provide more constraints on the EOS, which are 
correlated with various astrophysical observables. The 
FSUGold EOS is considerably softer than NL3 at high 
density, but stiffer at subnuclear density. The original 
FSUGold EOS supports a maximum NS mass of 1.7 M0, 
which is in contradiction to the recent observation of 2 
solar mass NS [2^ . We introduce a modification term to 
the original FSUGold interaction that increases the pres- 
sure at high densities so that the maximum NS mass is 
increased to 2.1 Mq. 

The paper is organized as follows. In section 2 we 
present the RMF models for NL3 and FSUGold effective 
interactions. Then a modification to the FSUGold effec- 
tive interaction is presented, which results in a maximum 
NS mass of 2.1 Mq. We also summarize numerical details 
for the interpolation scheme we use to generate the full 
EOS table. Section 3 shows results for our EOS, includ- 
ing various thermodynamic properties and the composi- 
tion. The r = beta equilibrium EOS is also presented 
and used to calculate the NS mass-radius relation. We 
compare our EOS and existing L-S and S-S EOS tables. 
Section 4 presents a summary of our results and gives 
an outlook for future work. In the appendix the public 
access to our EOS table is given. 



II. FORMALISM 



In this section we describe our relativistic mean field 
interaction, computational method, and interpolation 
scheme. 



Relativistic mean field interactions: FSUGold 
and NL3 



In this subsection we discuss the RMF interactions 
FSUGold and NL3. At low density the Virial expan- 
sion is the same in both EOS and has been discussed 
in a previous work [isj . We note there are more low 
temperature (< 1 MeV) data in Virial gas in the FSU 
EOS, which only induces small difference between NL3 
and FSU EOS. The comparison between NL3 and FSUG- 
old effective interactions has been discussed extensively 
in Ref. [2^. For completeness, we briefly present the 
RMF models to motivate this work. 

The RMF theory that we use is based on Ref. [Ill 
supplemented with an isoscalar-isovector coupling as in- 
troduced in Ref. 



25|. The basic ansatz of the RMF 



theory is a Lagrangian density where nucleons interact 
via the exchange of sigma- ((/)), omega- {V^), and rho- 
(bp) mesons, and also photons (v4^), given by [2^ [25t 
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(1) 



The values for various parameters in the interacting la- 
grangian can be found in Table I of Ref. [2^. Here we 
duplicate the bulk properties of infinite nuclear matter in 
Table HI The incompressibility, symmetry energy and its 
slope at nuclear saturation density predicted by FSUG- 
old are considerably smaller than NL3, i.e., the FSUGold 
EOS is much softer than NL3 at saturation density and 
above. Therefore FSUGold gives a smaller maximum NS 
mass, which will be discussed in following sections. How- 
ever, we note that FSUGold gives slightly stiffer nuclear 
EOS at subnuclear density than NL3. As a result, FSUG- 
old predicts a smaller neutron skin thickness in ^'^^Pb, 
0.21 fm compared to 0.28 fm given by NL3 The 
experimental constraints on these nuclear properties at 
best can only be inferred indirectly. The incompressibil- 
ity is constrained by giant resonance of nuclei, while the 
symmetry energy and its derivative can be constrained 
for example by studying the neutron skin thickness of 
^''^Pb, or comparing with microscopic calculations us- 
ing chiral perturbation theory or quantum Monte Carlo 
method. Both constraints on FSUGold and NL3 have 
been studied extensively in Ref. [23] . Again, FSUGold is 
consistent with these constraints. 

To better describe low density neutron matter, which 
is close to a unitary gas, we also introduced a density de- 
pendent coupling gs between scalar mesons and nucleons 
for FSUGold, which is the same as that for NL3 in Ref. 
p^ . Then we solve the Hartree problem at finite tem- 
perature for non-uniform nuclear matter at intermediate 
densities, and uniform nuclear matter at high densities, 
following Ref. [H. 
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Model 


po (fm-^) 


eo (MeV) 


Ko (MeV) 


J (MeV) 


L (MeV) 


NL3 
FSU 


0.148 
0.148 


-16.24 
-16.30 


271.5 
230.0 


37.29 
32.59 


118.2 
60.5 



TABLE I: Bulk parameters of infinite nuclear matter at sat- 
uration density Pg. The quantities and Ko represent the 
binding energy per nucleon and incompressibility coefficient 
of symmetric nuclear matter, whereas J and L represent the 
symmetry energy and derivative of the symmetry energy with 
respect to density at p^. 



B. Modified FSUGold effective interaction: FSU2.1 

The recent observation of a 1.97 ± 0.04 Mq NS [Hi 
puts a new constraint on the maximum mass allowed 
by effective interactions. The original FSUGold inter- 
action predicts a maximum NS mass of only 1.7 Mq. 
The FSUgold interaction is fitted to many observables 
near nuclear density, and the extrapolation to high den- 
sities is incompletely constrained. We consider a modifi- 
cation at high densities that increases the maximum NS 
mass. We assume additional repulsion between nucleons 
at short distances that is primarily isoscalar. This will 
lead to a correction to the pressure that is approximately 
independent of proton fraction. Furthermore, matter is 
degenerate at high densities so that the correction should 
be nearly independent of temperature. For simplicity we 
add the following term to the pressure. 



dP 



(2) 



where p is the baryon density and pq is a constant cor- 
responding to 0.2 fm"-^. Note that the extra pressure is 
zero for p < po- The corresponding change to the en- 
ergy is determined from the density integral of dP. The 
constant a=2 x 10^^ MeV^^ was chosen to increase the 
maximum NS mass to about 2.1 Mq. We refer to the 
EOS built with the original FSUGold effective interac- 
tion as FSU1.7, while the EOS with the modified effective 
interaction will be called FSU2.1. 

This procedure makes no corrections at low densities. 
This leaves unchanged all of the FSUgold predictions for 
nuclear observables. An alternative, and more consis- 
tent, procedure is to refit the effective interaction with 
a smaller ( coupling in Eq. (1). This will increase the 
pressure at high densities and lead to a larger maximum 
NS mass, l23j. An equation of state with this modified 
interaction will be presented in future work. However, 
we expect our simple prescription to give similar results. 



C. Computational methods 

We calculate the FSU EOS for many points in the pa- 
rameter spaces of density hb, temperature T, and pro- 
ton fraction Yp using the procedure in our previous pa- 
per [13]. There are a total of 107,000 points in the 3- 
dimensional parameter spaces of T, 1^, and ub- The 



overall FSU EOS took about 200,000 CPU hours and 
was run primarily on the Teragrid supercomputer cluster 
Ranger. Details about these parallel computations have 
been discussed in previous work [l3, [ij] . 



D. Bicubic interpolation 

As discussed in our previous paper [i3|, we use a hy- 
brid method of interpolation to generate an EOS table 
on a finer grid of density and temperature points. The 
thermodynamic pressure Pth can be obtained numerically 
from the free energy per baryon F/A, 
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For the zero temperature EOS, we smooth the pressure 
as a function of density for the calculated data points in 
Tableim Note the upper limit in densty is 10°-^ and 10°-^ 
fm^'^ for FSU2.1 and FSU1.7, respectively. This smooth- 
ing involves comparing the pressure at a given density to 
the pressure interpolated from neighboring points. We 
smooth the pressure in a 10 points per decade table by 
removing points that differ significantly (> 10"'^) from 
the geometric mean of neighboring points and replace 
them with interpolated one. Then we interpolate these 
smoothed pressures to create a finer table with 40 points 
per decade of density. The resulting energy at zero tem- 
perature is easily obtained by integration of these pres- 
sures with respect to density. The pasta phase has very 
close thermodynamic properties as normal nuclei [llj. 
Therefore we do not expect the smoothing procedure will 
lose important features of pasta phase. The purpose of 
smoothing is just to smooth the fluctuations in numerical 
evaluation of pressure and entropy. 



TABLE II: Range of temperature T, baryon density ub and 
proton fraction Yp in the coarse EOS table. 



Parameter 


minimum 


maximum 


number of points 


T [MeV] 


0, 10-"'' 




36 


logio(nfl) [fm^^] 


-8.0 


0.2 (0.4) 


83 (85) 


Yp 


0, 0.05 


0.56 


53 



Next, we calculate the entropy at finite temperatures 
T < 12.5 MeV. From the free energy at temperature T, 
density n b , and proton fraction Yp as in Table the 
entropy per baryon sth can be obtained numerically. 



Sth 



Finally, the energy per baryon eth is 

F/A - Tsth. 



eth 



(4) 



(5) 



For the finite temperature EOS with T < 12.5 MeV, we 
smooth the entropy, as a function of density, while en- 
suring thermodynamic stability. Then we perform mono- 
tonic cubic Hermite interpolation [2^ on this small table 
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of entropy with ten points per decade, to generate a larger 
table with 40 points per decade in both the temperature 
and density axes, as in Table lllll Now we integrate this 
smoothed entropy table as a function of temperature to 
get values of the free energy (adding the energy at zero 
temperature to make the full free energy). Thus we ob- 
tain the free energy on this finer 40 by 40 points per 
decade grid in a thermodynamically consistent fashion. 

Note that for higher temperatures T > 12.5 MeV, mat- 
ter is uniform for all proton fractions and densities. For 
uniform matter, all thermodynamic quantities can be ob- 
tained directly from RMF calculations, with good ther- 
modynamic consistency. 



TABLE III: Range of temperature T, density ub, and proton 
fraction Yp in the finely spaced interpolated EOS table. 



Parameter 


minimum maximum 


number of points 


T [MeV] 


0, 10'"'* lQ^>i'^ 


109 


logio(ns) [fm~^] 


-8.0 0.175 (0.375) 


328 (336) 


Yp 


0, 0.05 0.56 


l(Yp=0)-F52 



Finally, we carry out bicubic interpolation of the pre- 
vious free energy values (as in Table IIII[) to generate 
the entropy and pressure by thermodynamic derivatives 
Eqs. pi4p . This prescription guarantees the monotonic- 
ity of entropy and pressure in the final table, and con- 
serves the first law of thermodynamics in adiabatic com- 
pression tests. 

We first apply bicubic interpolation [g^l for the free 
energy. The first derivatives on the grid points are gen- 
erated from monotonic cubic Hermite interpolation [26|. 
The second derivative - the cross derivative d'^F/dndT, 
on the grid points is generated as in Ref. [131 . The bicubic 
interpolation can then fit free energies with cubic func- 
tions in temperature and density coordinates, and pro- 
vides the first and second (cross) derivatives. Then the 
entropy and pressure are obtained from Eqs. (|3l4p . Fi- 
nally, the boundary points at the highest temperature or 
density are discarded to avoid boundary artifacts in the 
interpolation. We note that our interpolation procedure 
is general and should work independent of the approxi- 
mations used to calculate the original free energy values. 
We test the thermodynamic consistency of our table in 

Sec. Impl 

The neutron chemical potential is derived from /i„ ~ 
d{nBF/A)/dnn\T,np, while /Xp = {F/A + p/ns - fJ.n{l - 
Yp))/Yj,. 



III. RESULTS 

In this section we present results for the free energy, 
composition, pressure, NS structure, and adiabatic index 
predicted by our FSU EOS. We compare our results for 
the FSU EOS together with NL3 EOS to those for the 
Lattimer-Swesty or H. Shen et al. EOS. 

A. Free energy and Composition 

In this subsection we discuss the free energy, average 
mass number of heavy nuclei, and mass fractions of dif- 
ferent components (n, p, a and nuclei) for the FSU EOS. 
In Fig. [1] the free energy per baryon F/A for matter at 
T = 1, 3.16, 6.31, and 10 MeV with different proton frac- 
tions is shown as a function of density ub- The free en- 
ergy F/A is obtained from Virial gas, nonuniform Hartree 
mean field, and uniform matter calculations. This figure 
resembles Figure 1 in our previous paper for the NL3 
EOS fl^. The transitions between the Virial expansion 
and the Hartree calculations (vertical dashed lines), or 
between the Hartree calculations and the uniform mat- 
ter calculations (vertical solid lines) are at the densities 
where Hartree or uniform matter calculations give the 
lower free energy. The transition between Virial EOS 
and RMF EOS is pretty smooth, as discussed in Ref. [l3|. 
However, there are errors in calculating the pressure and 
entropy numerically. It is the purpose of smoothing pro- 
cedure outlined in Sec. Ill Dl to remove these numerical 
errors. 

Figure [5] shows the average mass number A of heavy 
nuclei (with A > A) versus baryon density at different 
temperatures and proton fractions. Note that a parti- 
cles are not counted as heavy nuclei. Let us first look at 
the upper left panel (a), where T = 1 MeV. Nuclear shell 
effects give rise to several approximate plateaus in A vs 
density for each Yp, for example A ~ 50, ~ 80 and ~ 110. 
Usually A is larger in matter with smaller Yp. There are 
oscillations in A in the Hartree mean field regime. This 
is due to both nuclear shell effects and small errors in 
the free energy minimization due to our using a finite 
step in the Wigner Seitz cell size [l^ . The average mass 
A can be as large as 4,000 at high density. These large 
A values represent shell states or spherical pasta config- 
urations |11| . At higher temperatures, as shown in the 
other panels (b-d) of Fig. [2l A grows rapidly to several 
thousand in a narrow range of density before A drops 
abruptly at the transition to uniform matter. 



In Fig. [21 the mass fractions of unbound neutrons, un- 
bound protons, a particles and heavy nuclei are shown 
versus baryon density for a proton fraction of 1^ = 0.05. 
The solid lines are for the FSU EOS while the dashed 
lines are for NL3 EOS. Note that for Hartree mean field 
calculations there is only a single nucleus associated with 



each Wigner-Seitz cell. We define a nuclear level to be 
unbound when it has positive energy. The upper panels 
(a), (b) are for T = 1 and 3.16 MeV. The difference be- 
tween NL3 and FSU EOS appears around 3x 10"'' fm'^ 
for T = 1 MeV and lO"^ f^^-a f^^. y ^ 3 ^6 MeV, when 
the Virial gas to Hartree mean field transition occurs. 
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FIG. 1: (Color online) Free energy per baryon at temperatures of T = 1 (a), 3.16 (b), 6.31 (c), and 10 (d) MeV. The proton 
fraction ranges from Yp = 0.05 to 0.5. The transitions between the Virial expansion and the Hartree calculations are vertical 
dashed lines. The transitions between Hartree calculations and the uniform matter calculations are vertical solid lines. 



One can see that FSU and NL3 EOS differ in the mass 
fractions of heavy nuclei in the Hartree calculations. For 
T = 3.16 MeV, FSU and NL3 EOS also differ in the mass 
fractions of protons and a particles, due to different tran- 
sition densities. The bottom panels (c), (d) give the mass 
fractions of different species at higher temperatures, T = 
6.31 and 10 MeV. The differences between FSU and NL3 



EOS are small at high temperatures. Fig. 2] shows the 
mass fractions as were shown in Fig. [3J but for matter 
with higher proton fraction Yp = 0.5. Again, there are 
small differences between the FSU and NL3 EOS for the 
mass fractions, particularly for a particles, due to their 
different transition densities from Virial to Hartree. 



B. Neutron Star Structure 

Here we present the zero temperature EOS of nuclear 
matter in beta equilibrium and the resulting NS mass ra- 
dius relation. The results for the FSU1.7, FSU2.1, and 
NL3 EOS are compared. Fig. [5] shows the pressure versus 
density for zero temperature nuclear matter in beta equi- 
librium. The FSU1.7 and FSU2.1 EOS are stiffer (higher 
pressure) than the NL3 EOS for intermediate densities 



above 10~^ fm~^, and then softer than NL3 when the 
density is above 5xl0~^ fm~^. FSU2.1 becomes stiffer 
than FSU1.7 when the density is above 0.2 fm~^. The 
stiffness of the EOS at high density determines the max- 
imum NS mass. Therefore FSU1.7, FSU2.1, and NL3 
predict maximum NS masses of 1.7, 2.1, and 2.8 Mq re- 
spectively, as shown in the left panel (a) of Fig. [6l The 
FSU1.7, FSU2.1, and NL3 EOS predict radii of roughly 
13 km, 14 km, and 15 km, respectively, for a 1.4 Mq 
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FIG. 2: (Color online) Average mass number of heavy nuclei at temperatures of T = 1 (a), 3.16 (b), 6.31 (c), and 10 (d) MeV. 
The proton fraction ranges from Yp = 0.1 to 0.5. 



NS. The right panel (b) of Fig. [S] shows the NS mass vs 
central density. For a NS with given mass, the central 
density of the NS is inversely related to the stiffness of 
the EOS at high density. 



C. Adiabatic index 

In this subsection we discuss the adiabatic index F of 
nuclear matter, at constant entropy, for several differ- 
ent EOS. Figure [7] shows the temperature of the adiabat, 
with entropy S/A = 1, for nuclear matter with proton 
fractions of = 0.1, 0.2, 0.3, and 0.4. Note that baryon, 
electron, positron, and photon contributions to S are in- 
cluded. The sohd lines are for the FSU1.7 EOS while 
the dashed lines are for the NL3 EOS. Above saturation 
density FSU1.7 gives a slightly higher temperature than 
NL3. Since the modification term, Eq. ©, in FSU2.1 
does not depend on temperature, there is no difference 
between FSU1.7 and FSU2.1 in the adiabats. 



The adiabatic index F, 

describes the stiffness of the EOS at constant entropy. 
In Fig. m the adiabatic index is shown versus density 
for nuclear matter with constant entropy S/A = 1 and 
Yp 0.1, 0.2, 0.3, and 0.4. At subnuclear density, F 
has only small fluctuations with density. It then rises 
rapidly at the transition from nonuniform to uniform 
matter. This characterizes the stiffening of the EOS due 
to the large nuclear incompressibility. Beyond satura- 
tion density, NL3 has a larger value of F than FSU1.7. 
For FSUGold EOS, we incorporated more data for Virial 
EOS at low temperatures, which leads to small variations 
in the numerical evaluation of pressure and entropy at 
low density. This induces small differences between NL3 
and FSU EOS in the adiabatic index and temperature at 
low densities. 

In Fig. [S] we present the adiabatic index of various EOS 
assuming an initial entropy so that T = 0.7 MeV at a 
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FIG. 3: (Color online) Mass fractions of unbound neutrons, unbound protons, q particles and heavy nuclei versus baryon 
density at four different temperatures T = 1 (a), 3.16 (b), 6.31 (c), and 10 (d) MeV, and fixed proton fraction Yp = 0.05. The 
solid lines are for FSU(1.7) EOS while the dashed lines are for NL3 EOS. 




FIG. 4: (Color online) Mass fractions of unbound neutrons, unbound protons, a particles and heavy nuclei versus baryon 
density at four different temperatures T = 1 (a), 3.16 (b), 6.31 (c), and 10 (d) MeV, and fixed proton fraction Yp — 0.5. The 
solid lines are for FSU(1.7) EOS while the dashed lines are for NL3 EOS. 
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FIG. 6: (Color online) Mass versus radius, left (a), and mass versus central density, right (b), for NSs using FSU1.7 (solid), 
FSU2.1 (dot-dashed) and NL3 (dashed) EOS as shown in Fig. [1 
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FIG. 7: (Color online) Temperature of adiabat with S = \ for different proton fractions: (a) 0.1, (b) 0.2, (c) 0.3, and (d) 0.' 
Dashed (red) lines are for NL3 based EOS and solid (black) lines are for FSUGold EOS in this work. 



density ub = 10"'' fm-^ for Fp = 0.3. All EOS show a 
rapid rise of F beyond the transition to uniform matter ^ 
0.1 fm^'^, due to the incompressibility of nuclear matter. 
There are different versions of the Lattimer-Swesty EOS 
depending on the assumed value for the incompressibility 
of nuclear matter. We refer to LS180, LS220, and LS375 
for the Lattimer-Swesty EOS with Kq = 180, 220, and 
375 MeV, respectively. The LS180, LS220, LS375 and 
H. Shen EOS have a significant dip in F at densities just 
below the transition to uniform matter. In contrast the 
NL3, FSU1.7, and FSU2.1, EOS, do not show an obvious 
dip. 

The Lattimer-Swesty EOS is based on a simple liquid 
drop model and assumes a first order phase transition be- 
tween a low density vapor phase and a high density liquid 
phase. During this transition the pressure is indepen- 
dent of density so F is small. Indeed a first order phase 
transition may be appropriate for small systems such as 
heavy ion collisions. For small systems. Coulomb effects 
are relatively mild so that one can form a single uniform 
liquid phase. However, for large astrophysical systems. 
Coulomb interactions play a crucial role and prevent the 
formation of a uniform liquid phase. Instead, one must 
form a nonuniform phase with average proton density 
equal to the electron density. For this nonuniform phase, 
our relativistic mean field calculations give an adiabatic 
index that increases with density and does not show a 



decrease. Therefore, we believe that the decrease in F 
for the LS EOS may be an artifact of the approximations 
that they use and their assumption of a first order phase 
transition. The behavior of F for the H. Shen EOS may 
be related to their Thomas-Fermi approximation. For 
example, it is not clear if they considered a full range of 
possible shapes for nonuniform matter such as the shell 
states 11] that could be present in this density range. We 
emphasize that the H. Shen EOS involves variational cal- 
culations and Thomas Fermi approximations to the rel- 
ativistic mean field formalism that we calculate directly. 

The FSU1.7 and FSU2.1 EOS have a larger F com- 
pared to NL3 in the density range ub = 0.07 to 0.1 
fm~'^. This is probably related to differences in the den- 
sity dependence of the symmetry energy, see L in Table 
m Since L is smaller for the FSU EOS, it actually has 
a larger symmetry energy and pressure at low densities 
than NL3. This could lead to the larger F at low densi- 
ties. Nuclear measurements can constrain L and reduce 
the model uncertainty in F at low densities. For example, 
L can be constrained by measuring the neutron radius of 
^"^Pb with parity violating electron scattering j2^, [2^ or 
with measurements of isospin diffusion in heavy ion colli- 
sions [lO] • The possible impact of the density dependence 
of the symmetry energy on the EOS should be studied 
further. 

The rising of F for FSU2.1 at 71^=0.2 fm^^, due to the 
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FIG. 8: (Color online) Adiabatic index F along adiabat for different proton fractions Yp-. (a) 0.1, (b) 0.2, (c) 0.3, and (d) 0.4. 
Dashed (red) lines are for NL3 based EOS and solid (black) lines are for FSU1.7 EOS in this work. 



modification term to the pressure (and energy) beyond 
this density, indicates stiffening of EOS at high density. 
In a more consistent way to increase the pressure at high 
densities by refitting the effective interaction as discussed 
in Sec. IIB, F may rise more gently with density, which 
we will explore in future work. However the rise of F 
in FSU2.1 does not violate any physical principle we are 
aware of, especially the causality. Moreover, the peak at 
0.2 fm~'^ becomes smaller and the small kink disappears, 
if the EOS table has a lower resolution, for example 20 
points per decade in temperature and density. 

D. Adiabatic Compression Tests 

It is important that an EOS table is thermodynami- 
cally consistent. Otherwise small errors in the table can 
lead to artificial increases in entropy during an astro- 
physical simulation. Here we present adiabatic compres- 
sion test results for FSU1.7 and FSU2.1 EOS. The de- 
tails of this test have been discussed in a previous paper 
[l^ . The entropy should be conserved in this adiabatic 
compression independent of the density. In Fig. [TUl the 
temperature and entropy versus density during adiabatic 
compression are shown for nuclear matter with fixed pro- 
ton fraction Yp = 0.15 (left) and 0.50 (right). The initial 
density is 5.16738 x lO"*' fm'^ and initial T = 0.5 MeV. 



The FSU and NL3 EOS conserve entropy within 1%, ex- 
cept for a slight rise in S/A at extremely high density 
when yp=0.5. For a comparison test results for the H. 
Shen EOS is also shown (as red dashed curve) for sim- 
ilar initial conditions. Note that it is important to use 
an accurate interpolation scheme with the EOS table to 
ensure that the first law is satisfied so that entropy is 
conserved. The adiabatic compression test result for the 
H. Shen EOS was obtained using the routine developed 
in Ref. [3l|. In typical adiabatic compression tests, one 
has to interpolate the existing tables. For nuclear mat- 
ter at high densities with high Yp, the thermodynamic 
properties change fast as density rises. The NL3, FSU, 
HShen EOS used for this plot have 40 points per decade 
in density and temperature axes. For a finer table with 
80 points per decade, we find the entropy is conserved 
better at high densities (within 1 % deviation). However, 
this would increase the volume of table by a factor of 
four. Since this region is most likely not relevant in typ- 
ical supernova simulations, we decided to generate the 
smaller 40 points per decade table. 

E. Comparisons with existing EOS 

It is useful to compare the thermodynamics of the new 
FSU EOS with the NL3 EOS, the Lattimer-Swesty EOS 
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FIG. 9: (Color online) Adiabatic index F for different equations of state (EOS) assuming an initial entropy so that T — 0.7 
MeV at ns = 10"* fm"^. The proton fraction is Yp = 0.3. The FSU1.7 and FSU2.1 EOS are from this work and are identical 
below ub = 0.2 fm~'^, while NL3 is from ref. [ij and S-S is from ref. For the Lattimer-Swesty EOS [3] three curves are 
shown for different versions with incompressibilities Kq of 180, 220, and 375 MeV as indicated. 




FIG. 10: (Color online) Left panel: Temperature (a) and entropy (b) versus density in the adiabatic compression for nuclear 
matter with fixed proton fraction 0.15, for HShen EOS, NL3 EOS and FSU1.7 EOS. Right panel: same as left but with proton 
fraction 0.5. 
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FIG. 11: (Color online) Pressure versus baryon density from our EOS, Lattimer-Swesty's (L-S) and H. Shen et al.'s (S-S) EOS, 
with r = 1 MeV, Yp = 0.05 (a), T = 1 MeV, Yp = 0.4 (b), T = 6.31 MeV, Yp = 0.05 (c), and T = 6.31 MeV, Yp = 0.4 (d). 




FIG. 12: (Color online) mass fraction of heavy nuclei versus baryon density from our EOS, Lattimer-Swesty's and H. Shen et 
al.'s EOS, with T = 1 MeV, Yp = 0.05 (a), T = 1 MeV, Yp = 0.4 (b), T = 6.31 MeV, Yp = 0.05 (c), and T = 6.31 MeV, Yp = 
0.4 (d). 
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FIG. 13: (Color online) Average mass number A (black upper curves) and atomic number Z (red lower curves) of heavy nuclei 
for the FSU EOS (double dashed-dot), NL3 EOS (solid) , Lattimer-Swesty's (dot-dashed) and H. Shen et al.'s EOS (dashed) 
for r = 1 MeV, Yp = 0.4. 



(L-S), that uses a simple liquid drop model, and the H. 
Shen et al. EOS (S-S), that uses the Thomas Fermi ap- 
proximation and variational calculations with approxi- 
mating shapes for nuclei. The L-S EOS quoted in this 
subsection corresponds to the one with incompressibility 
coefficient i^o^lSO MeV, i.e. LS180. The comparisons 
for pressure, entropy, and heavy nuclei fractions between 
NL3, LS180, and S-S EOS have been discussed in a pre- 
vious paper p^ . In this subsection, we will focus on 
the differences between the FSU1.7 and the other EOS. 
The only difference between FSU1.7 and FSU2.1 EOS is 
the pressure and energy at densities above 0.2 fm^"^, as 
explained in Sec. Ill Bl 

In Fig. [11] the pressure for matter at T = 1, or 6.31 
MeV and Yp = 0.05 or 0.4 is shown for the FSU EOS, 
NL3 EOS, LS180 and S-S EOS. The pressure includes 
contributions from electrons, positrons, and photons. For 
neutron rich matter at subnuclear density as in the left 
upper panel (a) for T^l Mev and yp=0.05, the FSU 
EOS gives slightly larger pressure than all the other EOS. 
As the temperature and/or proton fraction grows, this 
difference becomes less obvious. On the other hand, for 



nuclear matter above saturation density, the FSU EOS 
gives the smallest pressure among all of these EOS. 

Fig-ttHcompares the mass fraction of heavy nuclei from 
the FSU EOS, NL3 EOS, LS180 EOS and S-S EOS, for 
matter at T = 1 or 6.31 MeV, and Yp = 0.05 or 0.4. 
The mass fraction of heavy nuclei can be very different 
among the various EOS and this causes the differences in 
entropy. However, the FSU and NL3 EOS only differ in 
the mass fraction of heavy nuclei for neutron rich matter 
at subnuclear density and low temperature as shown in 
the left upper panel for T=l Mev and yp=0.05. This 
is also discussed in Fig. [31 Note possible reasons for the 
differences in the mass fractions between NL3, LS180 and 
S-S EOS have been discussed in the previous paper [l^ . 
The sharp rise in the mass fraction for FSU and NL3 
EOS with T = 6.31 MeV and Yp = 0.4, also T = 1 MeV 
and Yp — 0.05, is due to transition from Virial to RMF 
results. For RMF, we solved the nuclear energy levels in 
non-uniform nuclei exactly and denoted nucleon as free 
when its energy is positive. This is different from Virial 
EOS and it is possible the composition is not smooth 
everywhere. The average A and Z of heavy nuclei for T 
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FIG. 14: (Color online) Entropy per baryon versus baryon density from our EOS, Lattimer-Swesty's and H. Shen et al.'s EOS, 
with T = 1 MeV, Yp = 0.05 (a), T = 1 MeV, Yp = 0.4 (b), T = 6.31 MeV, Yp = 0.05 (c), and T = 6.31 MeV, Yp = 0.4 (d). 



= 1 MeV and Yp = 0.4 is shown in Fig. [T31 This figure 
shows the effects of shell structure, that is included in 
FSU and NL3 EOS but is neglected in both LS180 and 
S-S EOS. As a result, A and Z for our EOS have a series 
of steps while A and Z for LS180 and S-S EOS increase 
smoothly with density. Note the NL3 EOS gives larger 
A and Z in the Hartree non-uniform region compared to 
the FSU EOS. NL3 also shows somewhat of a plateau 
in A and Z, where nuclei are in the shell shape [ll| . 
while the FSU EOS does not give shell shape nuclei for 
this value of Yp (it does for most values of 1^). More 
detailed discussion about shell shape nuclei can be found 
in Ref. 

Fig. [11] compares the entropy from FSU EOS, NL3 
EOS, LS180 EOS and S-S EOS, for nuclear matter at 
T = 1 or 6.31 MeV, and Yp = 0.05 or 0.4. Note that 
heavy nuclei have low entropy at low temperatures so 
the difference in entropy between the FSU and NL3 EOS 
is small. While at higher temperatures, the entropy is 
dominated by unbound nucleons and photons and the 
two EOS agree very well. 

IV. SUMMARY 

In this paper we present a second equation of state 
(EOS) of nuclear matter for a wide range of tempera- 
tures, densities, and proton fractions for use in super- 
novae, neutron star - neutron star mergers, neutron star - 



black hole mergers, and black hole formation simulations. 
We use a density dependent relativistic mean field (RMF) 
model, FSUGold, for nuclear matter at intermediate and 
high density with a spherical Wigner-Seitz approxima- 
tion for nonuniform matter, which incorporates nuclear 
shell effects [iSj. For nuclear matter at low density [is!], 
we use a Virial expansion for a non-ideal gas consisting 
of neutrons, protons, a particles and thousands of heavy 
nuclei from the finite range droplet model (FRDM) mass 
table [32J . The difference from the first EOS is that here 
we use the RMF effective interaction FSUGold whereas 
our first EOS was based on the NL3 effective interaction. 
The Virial gas part is common and has been discussed in 
a previous paper (Tsj . 

We tabulate the resulting EOS at over 100,000 grid 
points in the temperature range T = to 80 MeV, the 
density range ub = 10~^ to 1.6 fm~^, and the proton 
fraction range = to 0.56. We present differences be- 
tween our FSUGold based EOS and our NL3 based EOS 
along with some existing EOS, for the thermodynamic 
properties, composition, and neutron star structure. In 
particular, we focus on the differences between the FSU 
and NL3 EOS. The FSUGold EOS is considerably softer 
than NL3 EOS in the high density, but stiffer at subnu- 
clear density, as shown in Fig. [TT] The FSU EOS also 
predicts different mass fractions and average A and Z of 
heavy nuclei from NL3 EOS, as shown in Figs. [31 [H and 

m 

As for the NL3 EOS, we use a hybrid interpolation 
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scheme to generate a full table for the FSU EOS on a fine 
grid that is thermodynamically consistent. This ensures 
that the first law of thermodynamics is satisfied and that 
entropy is conserved during adiabatic compression. 

Moreover, the original FSUGold EOS has a maximum 
neutron star mass of 1.7 Mq. We introduce a modifica- 
tion in the pressure at high density, which accommodates 
a 2.1 Mq neutron star as shown in Fig. [51 The EOS based 
on the original FSUGold is therefore called FSU1.7 while 
the modified EOS is called FSU2.1, according to their 
different maximum NS masses. Note in the final EOS ta- 
bles the upper limit in densty is 10°-^''^ and lO'^ '^'^^ fm^"^ 
for FSU2.1 and FSU1.7, respectively. 

Finally, the EOS tables for FSU1.7 and FSU2.1 are 
available for download as described in Appendix I VII Our 
goal in constructing EOS tables is two fold. First, we 
intend to calculate EOS tables with different pressures. 
This will allow one to correlate features of astrophysi- 
cal simulations with EOS properties. In this paper we 
calculate a new EOS that has a lower pressure, at high 
densities, than our previous NL3 EOS. Second, we aim to 
provide detailed composition information in future EOS 
tables. This can be important for neutrino interactions. 
Historically, most EOS tables, used in astrophysical sim- 
ulations, have only provided mass fractions for neutrons, 
protons, as, and a single average heavy nucleus. How- 
ever, electron capture on a range of heavy nuclei may 
be important for the proton fraction Yp during the infall 
phase of a supernova [34!|. In addition, mass 3 and other 
light nuclei can be important for antineutrino opacities 
[35I l36j . Deuterons have been included in current Virial 
EOS via the second Virial coefficient between neutrons 
and protons [l3l |. and we plan to include mass 3 and 
other light nuclei in future work. Note that the EOS ta- 
ble described in this paper only provides a single average 



heavy nucleus. This is done to be compatible with exist- 
ing simulations and neutrino opacity codes. However we 
will make available additional composition information 
in future work. 
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VI. APPENDIX: FORMAT OF EOS TABLES 

Here we describe where the tables 
FSU1.7EOSbl.01.dat FSU1.7EOS1.01.dat, 
FSU2.1EOS1.01.dat and FSU2.1EOS1.01.dat can 
be downloaded. The entries in these tables are the same 
as those in NL3 EOS tables and explained in a previous 
paper 0j. One should download the gzip compressed 
files (that are about 100 MB each) and use gunzip to 
decompress them. The grid structures of these tables 
are indicated in Table IIIII and contain approximately 
517 MB of data each. The tables, a sample FORTRAN 
computer program, and a readme file are available for 
download both at our website 



'http : / / Cecelia. physics . ind.iajia.edu/gajig_shen_eos/ 



and, from the Electronic Physics Auxiliary Publication The entries in the table are the same as that for the NL3 
Service - EPAPS web site Please check our web site EOS and have been discussed in the appendix of Ref. [lij . 
for any updated information regarding these EOS tables. 
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